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The lattice discrepancy of certain three-dimensional bodies 

Ekkehard Kratzel and Werner Georg Nowak M (Vienna) 



Abstract. Based on a very precise approximation to the lattice discrepancy of a Lame disc, an 
asymptotic formula is established for the number of lattice points in the three-dimensional body 



o 



o 

CO 



> 



o 
o 



|ui| mfe + (|u 2 | fe + |« 3 | fe ) m <* mfe , 



for large real a; and fixed reals m, fc . Particular attention is paid to the boundary points of Gaussian 
curvature zero. 
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1. Introduction. For a compact body B in R 3 the lattice point discrepancy of a 
copy xB, obtained by linear dilation by a large real parameter x, is of interest: 

P B (x) := # (xBnZ 3 ) - vol{B)x 3 . (1.1) 

The theory of its estimation, resp., asymptotic evaluation has been described in detail 



■<^j- ; in Kratzel' s books [13], [15], even in the frame of a more general s -dimensional setting, 
and also in Huxley's monograph [7] which concentrates on the planar case. See also the 
■ survey article [11] by the authors with Ivic and Kuhleitner. 

The case of a convex body C of smooth boundary dC with bounded nonzero Gaussian 



curvature throughout is quite well understood. Hlawka's [6] classic bound Pc(x) <C x 3 / 2 
has been refined up to 0(x 63 / 43+e ) by W. Miiller [23], using exponential sums in a most 
ingenious way. Furthermore, we also know that 



x 



ft W = n (*(log,)V), and yWdz«*3 (log ^. 



See Nowak [24], and Iosevich, Sawyer & Seeger [9]. For the special case of the sphere 
in M 3 , cf. Heath-Brown [5], for rational ellipsoids Chamizo, Cristobal & Ubis [1], for 
ellipsoids of rotation the present authors [19]. 

As soon as the boundary dB has points of curvature zero, the problem becomes much 
harder and, consequently, our knowledge is rather fragmentary. For partial results, see 
Haberland [4], Kratzel [14], [16], [17], [18] and Peter [26]. A fairly general theorem on 
bodies of rotation has been established by Nowak [25]. The method used there can 



The authors gratefully acknowledge support from the Austrian Science Fund (FWF) under project 
Nr. P20847-N18. 



2 



be called the "cut-into- slices approach": Since xB intersects every plane orthogonal 
to its axis of rotation in a circular disc, a very accurate approximation to the latter's 
discrepancy is employed, followed by a ("careful") summation with respect to the third 
coordinate. This idea has also been pursued by the authors for the special case of a 
rotating Lame's curve [20]. 

As a rule, with this sort of problems one cannot expect results which are both very 
sharp and general. Since the methods used are highly technical, in order to get precise 
estimates one can only deal with a particular class of bodies at a time. For the expert 
it will be transparent that the argument developed may be applied in other similar 
situations as well. 

On the basis of this understanding, in the present paper we will concentrate on bodies 

B m ,k := {(«i,«2,« 3 ) e K 3 : \ui\ mk + (\u 2 \ h + |« 3 |T < 1 } > ( L3 ) 
where k > 2 and to > 1 are fixed real numbers. 

Let A m ^(x) the number of integer points in xBm tk , then for integers k and to, 
A m ,k(x) describes the average number of representations of integers by the form \u\ | mfc -|- 
{\u 2 \ k + \u 3 \ k ) m : 

A m , k (T 1 /^)= 5> m , fc (n), 

n<T 

where 

r m ,k( n ) : = #{(ni,n 2 ,n 3 ) G Z 3 : \ ni \ mk + (\n 2 \ k + \n 3 \ k ) m = n} . 

The cases to = 1 (super sphere) and k = 2 (rotating Lame disc) have been dealt with 
by Kratzel [13], [14], resp., by the authors [20]. 

An important issue of the present paper is that the cut-into-slices method will be applied 
for the first time to a body which lacks symmetry of rotation. Thus, the planar discs 
arising are no longer circular, but they are bounded by Lame curves. The latter also 
contain points of curvature zero, which renders the situation much more complicated. 
We will establish sort of a "truncated Hardy's identity for Lame discs" which should be 
of some interest of its own - see Theorem 2 below. 

The body xB mjk has been investigated earlier by Kratzel [18], by a purely "three- 
dimensional" method: This lead to the estimation of double trigonometric sums and to 
a somewhat less precise final result. 

In fact, the Gaussian curvature of dB m ^ k vanishes on each curve of intersection with one 
of the coordinate planes. Where any two of these curves meet, there arise flat points, 
i.e., in (±1,0,0), (0, ±1,0), and (0, 0, ±1). We will be able to evaluate precisely the 
contribution of these flat points to the lattice discrepancy. The rest of the curves of 
Gaussian curvature zero will be taken into account by appropriate O -terms. (Cf. the 
remarks after the statement of the main theorem.) 

Theorem 1. For fixed reals k > 2, to > 1, with mk > |, and large real x, the 
number A mjk (x) of lattice points in the body xB m ^ satisfies the asymptotic formula 

A m ,k{x) = vol(B m ,k)x 3 + H mkjkjl (x) + H mkjkj2 (x) + K m , k (x) , 
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where 

H mW {x) = F mW {x) x 2 ~ 2 / mk + 0(x) , 

H mW {x) = T mW {x) x 2-Vmk-l/k + Q {x) ^ 

with continuous periodic functions J r m k,k,i(x) , i = 1,2. The functions H and T will 
be defined and analyzed in section 2. The remainder term is given by 

/ q 7 / 9 c\ _ / 339 131 , 18627mfc-20614 \ 

Um,k(x) = O [X 6 ' /Zb \ + O ix 2 08 104 m fc (\0gx) 8320T^fe j 

_/ 339 235 18627 (-> l / u\\ 

+ Olx 20S 208fe (log x) 8320 1 1 '). 

Thus, 

( x 37 / 25 if k < = 7.54 ..,andmfc<^ = 8.408 . . ., 

n ( \ ^ J 339 131 18627mfc-20614 ^cn , 2 R2 -i 11 

rCm,k\x) <€. < X208 i04 m fe(l O gx) 832 0m fe if mfc > and m > |§| = 1.11 . . ., 
^ xllK2w(log£)iMir (1 ~ 1/fc) if A; > and m < §§§. 



Remarks. As hinted at above, the term J^mk^^ix) x 2 2 / mk comes from the points 
(±1, 0, 0) where the two Lame curves 



|„, \mk i |„, \mk i \ |„, \mk i |„, \mk i 

£l . j Nl + m -i, and r „ . J Nl + NI -i 



intersect. The curvature of either of these curves vanishes at (±1, 0, 0) , of order mk — 2 
(related to the arclength). Similarly, the points (0, ±1,0), (0,0, ±1) which lie at the 
intersection of C\ , resp., Ci with 

N| fc + \u 3 \ k = 1, 
ui = 0, 

contribute the term Fmk,k,2(x) x 2 ~ l / mk ~ l / k . The curvature of £3 has zeros of order k — 
2 at (0, ±1, 0) , (0, 0, ±1) . As the proof will show, the contribution of the rest of C\ , £2 

to the lattice point discrepancy can be bounded by O (^x^~ wilk (logx) 86 832o m fe 06 j 
while £3 without the flat points contributes at most O (^x^~ ^ (logx)^^^ 1-1 ^^ . 



2. Auxiliary results. 

Lemma 1. (Vaaler's approximation of fractional parts by trigonometric polynomials.) 

For arbitrary w E M. and any integer H > 1 , let ip(w) := w — [w] — | , 

ip H (w) := - ^ a h,H sin(2nhw) , ip* H (w) \= ^ &h,H cos(2nhw) , 

0<h<H 0<h<H 
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where, for h = 1, . . . , H — 1 , 
an<i 

P (o = ^(i-ocot(7ro+^ (o < e < 1) 

TTien t/ie following inequality holds true: 

\ip{w) - V'hHI < i)* H (w) + ^ . 



Proof. This is one of the main results in Vaaler [27]. A very well readable exposition 
can also be found in the monograph by Graham and Kolesnik [3] . 

Lemma 2. (A) Let F e C A [A,B], and suppose that, for positive parameters 
X,Y,Z, we have 1 < B - A < X and 

< X 2 - j Y- x for j = 2,3,4, \F"\ > c^Y' 1 , 

throughout the interval [A, B] , with some constant Co > . Let J' denote the image of 
[A, B] under F' , and F* the inverse function of F' . Set further 

Jo ifF'^eZ, 

r(0 ■— jmin^lF'^ir 1 ,^) else, 

where \\-\\ stands for the distance from the nearest integer. Then, with e(w) = e 27 "™ as 
usual, 

£ e(F(n)) = e (^P) E * 1 e(F(F*(£)) - IF*(1)) + 



A<n<B 



+ 0(r(A)) + 0(r(B)) + O (log(2 + length( J'))) , 

with the notation 

a ^ n ^ b a<n<& 

where xz is the indicator function of the integers. 

(B) Let further G E C 2 [A,B], with < X~i Z for j = 0,1,2, then it follows 
that 



+ (z(Vy + log(2+ length(J')))) • 
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Proof. Transformation formulas of this kind are quite common, though often with 
worse error terms. For part (A), see Lemma 2.2 in Kiihleitner and Nowak [22]. Part (B) 
can be found as f. (8.47) in the recent monograph [12] of H. Iwaniec and E. Kowalski. 



The functions H a ^ t i(x) and H a ^ 2 {x). 

For what follows, compare Kratzel [18], and also his monograph [13]. We start with the 
definition of generalized Bessel functions 

l 



for reals rj > 1, v > 1/rj , and real x > . Further, let 

^J>\x) = 2yfrT[y + 1 - 1/ V ) (—) 4 v) (2™x) , 



n = l 



the series converging absolutely for v > 1/rj. On the basis of these functions, we define 
for reals a > b > 2 , 



and 



^ i(x) = Sxw^ (x) ' (2 - 1} 
i 

H aA2 (x) = 8x J t a -\l-t a ) l ^- l ^ h {xt)dt. (2.2) 
o 

Of course it is important to know the asymptotic behavior of these functions. Since 
an asymptotics for the generalized Bessel functions is provided by [13, Lemma 3.11], it 
readily follows that 

/ \ ^ / ,\ 2-2/0 sin(27rnx - ir/a) _ 1/ , 

H aA1 (x) = C 1 (a,b)x 2 2/a >Z n i+2/a + °( x )> 

n=1 " (2.3) 
2r 2 (l/6) 2 / a \2/« / 2^ 



To expand H a ^,2{x), we approximate J^^ymxt) again by [13, Lemma 3.11], deriving 
J<«(2™*) = -±= (^) '"in (a™* - + ((«*)-') . 

1/2 

at least for sufficiently large. The part f in the definition of H a ^^(x) is only an 



Oil) . The arising new integrals 



1 

J t a - 1 / 6 (l-t a ) 1 / a - 1 sin(2 7 rna;t- ^) dt 



1/2 
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are dealt with according to E.T. Copson [2, p. 24, formula (11.6)]. This finally yields 
the asymptotics 

oo 

W = C 2 (a, 6) x 2 " 1 /"- 1 / 6 ^ n- 1 " 1 /-!/^ sin ( 2 ™ - ^ - ^) + 0(x), 

n=l 



^ , ,v 16 a 1 /^ 1 / 6 / 1\ / 1 
C 2 (a,b) := — - — -n -ttt T 1 + - T 1 + - 



(2.4) 



3. Preparation of the estimate. The " cut-into-slices approach" gives 

A m ,k(x)= ^((^ mfc -i^r fc ) 1/m ) , (3.1) 

|m I <x 

where, for real W > , 

L fc (W) = X) 1 

|n 2 | fc + |n 3 | fe <W 

is the number of lattice points in a Lame disc with length parameter W x l k . According 
to [13, formulae (3.57) and (3.47)], 

L k (W) = a k W 2 ' k + 8I k (W) - 8A k (W) + 0(1) , (3.2) 

where a k = ^-p^/k) * s ^ e area °^ the urn ^ Lame disc, 

w 1/k 

h(W):= I ^)d((W-« fc ) 1/fc ) , 

' (3.3) 

A k (W):= Yl ^{{W-n k )^ k ). 

(^W) 1 / k <n<W 1 / k 

By Euler's formula, 

I k (W) = \ Yl ( W ~ \ n \ k Y /k -\ I (W -\u\ k ) 1/k du 
W k < w \u\J<w 



1 (W-\n\ k f/ k -^W 2 / k . 

\n\ k <W 



Combining this with (3.1) and (3.2), we obtain 

A m , k (x) = -a k (x mk -\ ni r k f imk) + 4S(x) 

|ni I <x 

-8 Y A ^ {^ mk - Ki mfc ) 1/m ) +o(x), 



\m I <x 
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with 

S(x):= ((a; m *-|n 3 r fc ) 1/m -|Ti 2 | fc ) 1/fc . (3.4) 

\n 2 \ mk + \n 3 \ mk <x mk 

According to [13, Lemma 3.12], 

|ni |<3! 



XI + A) ' 3/(mfe) (3.5) 

= vol(iB m> fc) x 3 + H mkjkjl (x) , 

with an appeal to the definition (2.1). Hence we arrive at 

A m ,k( x ) = -vol(B m , fe ) x 3 - H mkjkjl (x) + 4S(x) 



8 ^ A fc [{x mk - |ni| mfc ) 1/m ) +0(x). 



ni I <x 



(3.6) 



This formula already outlines the strategy for the completion of the proof: We will have 
to evaluate S(x) , and then to estimate the multiple fractional parts sum. 



4. Evaluation of the sum S(x) . 



Proposition. For large real x and fixed real numbers m > 1 and k > 2 , the sum 
S(x) defined in (3.3) satisfies the asymptotic formula 

S(x) = -vo\(B rn , k )x 3 + -H rnW (x) + -H mW (x) + ^ii-^(logx)™ (1 - 1/fc) ) . 



Proof. Obviously, 



\s{x)= E" Ux^-n^f^-ni) 1 '" = Efdh 

4 n mk +n mk^ x mk V 7 (t fc + „fe)r7i + „mfe < 



where E means that terms corresponding to n 2 n 3 = get weight ^ > an d the term 
with n 2 = n 3 = gets weight | . We subdivide this sum into 6 subsums, according to 
the relative size of the three variables of summation, resp., integration, t±, n 2 ,n 3 . Thus 
£1,2,3(20 comprehends the case t\ < n 2 < n 3 , and so on. Terms with n 2 = n 3 get 
weight I , which is symbolized by the notation E . Thus, 

5 1,2,3 = E / d*l 

{t k +n k ) m + n ™ k <x mk 
< t ^ <^2 ^ n 3 



= St 2 , 3 (x) + -P 2 M^) + 0(x), 
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where 



(.t k +t k ) m +n™ k <x mk 
0<t 1 <t 2 <n 3 



and -P 2)3 (:r) has been defined in [18, formula (9)], with (k, k) instead of (mk,k). It 
satisfies the estimate 



339 235 



>(18627/8320)(l-l/fc) 



P2,3( X ) ^ x2 ° 8 2 ° 8k Q°g x ) 

according to the first formula on top of page 768 of [18], if one replaces M. Huxley's 
bound in [7] by his more recent slight improvement in [8]. Further, by the Euler 
summation formula, 



Sl,2, 3 (*)= fff d(t lt t 2 ,t 3 ) ~ 



(t k +t k ) m +t™ k <x mk 
0<*1 <t 2 <*3 



J J V ((x mh - {t\ + t k 2 ) m ) x /( mfc )) d(ti, £ 2 ) + O(x) 

Q< t mfe< t mfe< xm fe_(' t fe_|_ t fe- )m 

+ ^H mW (x) +0(x), 



(t k +t k ) m +t™ k <x mk 
0<t 1 <t 2 <t 3 



by [18, formula (12)] and the formulae at the bottom of p. 762 of [18]. Thus, altogether, 

«Sl,2,3(aO= /// d(ti,t 2 ,t 3 ) + ^^rnfc,fc,l( ;r )+ (^ il "^( l0 g^)™ (1 " 1/fc) ) • 



(t k +t k ) m +t™ k <x mk 



0<ti<t 2 <*3 

(4.1) 

Similarly, 

Si, 3 ,2(x) = S* w (x) + ^PsA x ) + 0(x) , 

where 

Si,3,2(*)= E//d(*i,*3) 

(t k +n k ) m +t™ k <x mk 
0<t 1 <t 3 <n 2 

and ^3,2^) has been defined in [18, formula (16)], with (k,n) instead of (mk,k). It 
satisfies the same estimate as -P 2 ,3(rr) , appealing to [18, Lemma 2]. Consequently, again 
by Euler's formula, 

Si l3l2 (aO= /// d(fi,f 2 ,f 3 ) + Hl 2 (x) + 0{x^-^(l 0gx )W^(^/^ , 

(t k + t k ) m +t™ k <x mk 
0<*l< t 3< t 2 

(4.2) 
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with 

*;,(«) = - // ^(((^- ( ?')^-^)d (tl , t 3) 

0<t k <t k <(x mk -t™ k ) 1 / m -t k 

Moreover, 

S2,i,s(x) = SI 1>3 (x) + P^x) + O(x) , 

where 

5 2 V, 3 (x)= E//d(ti,t 2 ) 

(tj+t|) m +™^ fe <a: mfe 
0<*2 <*1 <™3 

and 

P 1 * 3 (x) = - E/ d ^i i> U(x mk -n™ k ) x l m -t\) x l k \ 

2- m (x mk -n™ k )<t™ k <x mk -n™ k 

The substitution 

z= {{x mk -nf k ) 1/m -t k ) 1/k & tl = {(x mk -nf k ) l/m -z k ) 1/k 

yields 

P U( X ) = £ / dz ^(z)z k - 1 Ux mk - n™ k y/ m - z k \ 1+1/k 

0<z mk <2- m (x mk -n™ k ) 

Hence, integrating by parts, we conclude that 

P* 3 (x) «x. 

Furthermore, essentially repeating an earlier argument, 



5 2 *i >3 (x)= JJJ d(t u t2,t 3 ) - 

(t k +t k ) m +t™ k <x mk 

0<*2 <*1 <*3 

^(i m *-(t{ + ft m ) 1 /W) d(*i,* 2 ) +0(x) 

0<t^ fe <t™ fe <x mfe -(tj+t^) m 

d(ti,t 2 ,t 3 ) + Y^H mkjk;l (x) +0(x) . 

(t k +t k ) m +t™ k <x mk 
0<t 2 <t 1 <t 3 




Hence also 



S$,i,3(x)= JJJ d{t u t 2 ,t 3 ) + L Hmktkjl ( x ) + 0(x). (4.3) 

(t k +t k )m +t rnk< x mk 
°<*2 <*1 <*3 
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Next, 
with 



and 



55,3,1^)= E//d(ti,t 2 ) 

0<t 2 <n^<t 1 



P£,i(*) = - E/dti v(((^ mfc -< fc ) 1/m -t? 

(• x mfe_ n mfe)l/ m _ n fe< t fe<(- x mfe_ n mfe- ) l/m 



Again by [18, Lemma 1], we bound P^i(x) by O [xzol ioik (log a?) 1 8 8 32o 7 ( 1 1 / k )j as well. 
Further, 

<S 2 * 3 ,i(*) = HI d(fi,f 2 ,f 3 ) + HsV^ + OCx), 

{tk+tk)m +t mk < x mk 
0<*2<*3<*1 

where 

#3*l(*) = - || V> ((^ " (t? +^) m ) 1/(mfc) ) dfa.fc) 

o<t™ fe <x mfe -(t* , +t^) m <t5 nfe 
|| ^(^2 mfc - 1 (x mfc -2 mfc )- 1+1/m ((x mfc -2 mfc ) 1/m -^)- 1+1/fc d(^,t 2 ), 

Q<t k <Z k <{x rnk -Z rnk ) 1 / m -t k 

on the basis of the substitution 

z = ( x mk -(t k 1 +t^ m ) i/{mk) & t 1 =((x mk -z mk ) i / m -t%y /k . 

Observing that t\ < ^(x mk — /2 mfc ) 1 / m and z mk < \x mk , an integration by parts shows 
that H^ t (x) <i. Hence, altogether, 

52,3,1(3:)= HI d(fi,f 2 ,f 3 ) + o(x«»-^(logx)W(i-i/*)) . (4.4) 

{t k +t k ) m +t™ k <x mk 
0<t 2 <*3 <t 1 

Further on, 

S3,i,2(x) = Sl h2 (x) + Pl 2 (x) + O(x) , 

with 

Sti,*(*)= E//d(ti,t 3 ) 

°<*3 <*1 <™2 
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and 

P* 2 ( x ) = - E / d*i i> (( xmk ~ (*i + ^) m ) 1/(mfc) ) • 

a! mfc_ t mfc<( t fe +n fc)m< !1 .mB' 
0<tj <«2 

In order to bound Pf 2 (#), we substitute 

z= (x mfc -(tf + n^) m ) 1/(mfc) ^ h = ({x mk -z mk Y' m -n 1 ^ 1 ^ 
and obtain 

P 1 * 2 (x) = E / d ^ ^(^)^n 2 (^) , 

(z k +n k ) m <x mk -z mk <2 m n 2 nk , z>0 

where 

^(s) = z mk -\x mk - z mfc )" 1+1/m ((z mfc - z ™ fc ) 1/m - n^~ 1+1/fc . 

The condition (z k + n 2 ) m < x mk — z mk on the one hand implies that z mk < \x mk , 
hence (x mfc - z ™ky\+i/m x ^fc-mfc _ Q n the other hand, (x mk - z ™ k fl™ > z k ^n k 2l 
thus ((x mk - z mk ) x l m - nl)~ 1+1 / k < z x ~ k . Therefore, altogether, 

<p n2 (z) <€.x k - mk z mk - k < 1, 

uniformly in n 2 . Since <fi n2 ( z ) increases monotonically with z, integration by parts 
readily gives 

P^(x) «x. 

Further, 

<S 3 * 1)2 (:r)= HI d(fi,f 2 ,f 3 ) + ^r, 2 (x) + 0(x), 

(t k +t k )m +t mk < x mk 
0<t 3 <t 1 <t 2 

where 

h; >2 (x) = - II v>(((* mfc -CY /m -*i) 1/fe ) d(ti,t 3 ). 

0<t™ fe <t™ fe <2~ m (x mfe -t™ fe ) 

Collecting results, we arrive at 

£3,1,2(3;)= HI d(t lt t 2 ,t 3 ) + HZ j2 (x) + 0(x). (4.5) 

(t k +t k ) m +t™ k <x mk 
0<t 3 <t 1 <t 2 

Finally, 

Sa,2,i(z) = SUJx) + P^rr) + 0(x) , 
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with 

S 3 *,2»= E//d(*l,*3) 

(t^+nk)m +t mk < x mk 
0<*3 <"2 <*1 

and 

tff,i(z) = " E / d*! V fe m * - (*} + nj)™) 1 /^) . 

x mfc_„mfe< {t fe + „fc )m < a , m fe 
0<"2 <*1 

Quite the same analysis as used before for Pi 2 0*0 applies again and yields 

P 2 *i0*0 «x. 



Furthermore, 



<S 3 * 2)1 0*0 = HI d(fi,f 2 ,f 3 ) + H^ + OCx), 



( t k +t k )m+t rn.k< x rnk 
0<t 3 <t 2 <*l 



where 



H 2 * = - || V (((^ mfc - #*) 1/m - if) 1 '*) d(ti,t 3 ) • 

0<t%<(x mk -t™ k ) 1 / m -t k <t k 

We transform this double integral by the substitution 

2 = - *™ fc ) 1/m ~ *i) ^ ti = ((x rofc -^*) 1/m -z*J 

to obtain 



# 2 * = || VOO-^" 1 ((^ mfc - #P*) 1/m - z k ^j 1+1/k d(U 



-1+1/fc 

0<^<Z fe <-§(x mfe -t™ fe ) 1 A 

Integration by parts shows again that H^^x) <C x, hence 




0*0= /// d(*i,* 2 ,* 3 ) + O(x) 



( t fc_|_ t fc)m_(_ t mfe ^^mfe 
0<i3<*2<*l 

Adding up the results (4.1) - (4.6), we get 

\s{x) = ^vol(B m , k )x 3 + ^H mW (x) + Hl 2 (x) + Hl 2 (x) 



_ / 339 235 18627 /■. -, /, \\ 

+ CM £ 2 08 208fe (log X) 8320 ^ ' ' . 
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Finally, 

HU*) + H Ux) = - ff ,(((^- t -)^-^)d( t M3). 

meix(2 m t™ k ,(t k +t k ) m )<x mk -t™ k , ti,t 3 >0 

According to [18], formulae (15) and (20), and p. 763-764, this last double integral 
equals -^H mk ^ k ^.{x) . This observation completes the proof of the Proposition. 

We conclude this section by the remark that using the Proposition in (3.6) gives 

A m ,k(x) = yol(B mj k)x 3 + H mk)ky i{x) + H mk)ky2 {x) 

-8 Afe^-Kr^^+O^ll-mQog^Sd-V^ . 



Ini I <x 



(4.7) 



5. Approximating the lattice discrepancy of the Lame disc. The lattice point 
discrepancy of the disc \u\\ k + \u2\ k < W is expressed by (3.2). In fact, for I k (W) a 
very precise evaluation is known, namely 

w = \ (s) "' r i 1 + 1) w^wEr 1 - 1 '- {^jw 1/k - + on) . 

See [13, p. 148]. For A k (W) as defined in (3.3), the sharpest upper bound known reads 

/ 131 18627 \ 

Afc(W) = O (w^(logW)-ssss-j . (5.1) 

This follows by the argument of Kuba [21], if one uses Huxley's method in its most 
recent form [8]. Our present target will be to approximate A k (W) by trigonometric 
polynomials, with a fairly small O-term. This result is intended to be useful for further 
calculations involving A k (W) , like a summation with respect to a third dimension, in 
the sense of (3.1). 

To this end, recall the notation of Lemma 1 and let H\y denote an arbitrary map from 
{n G Z + : t;W < n k < W } into the integers exceeding 1, so that max(logif^) ^ 
log IF. Further, define A k ^ Hw (W) and A* k Hw (W) analogously to A k {W) , replacing 
tp by ip Hw (n) , resp., ^* Hw(n) , for every value of n in (5.1). By Lemma 1, 

\A k (W) - A k , H (W)\ < A* kiH (W)+ J2 7T~Tn) ' ^ 

(±W) 1 / k <n<W 1 / k V ' 

Put for short a = (<x h ,H w (n)) 0<h<Hw(n y P = (Ph,H w (n)) 0<h<Hw{ny and write 7 = 
(lh,H w (n)) 0<h<Hw(n) for either a or /3. In fact, A kMw{n) (W) and A* k Hw{n) {W) can 
be transformed by the same calculation, considering exponential sums 

E, := E E lh,H w (n) e (-h(W - n*) 1 /*) . 

(%W) 1 / k <n<W 1 / k 0<h<H w (n) 
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We divide the range of n by a dyadic sequence N < Ni < . . . < Nj defined by the 
condition 



d 

du 



(-h{W -u k ) 1/k } 



= hV <^ N j = W 1/k (l + 2- j9 )- 1 ' k , (5.3) 

u=Ni 



with 

k 



k-1 

for short, throughout what follows. J is chosen as 

"log((W/c ) (1 - A) / fc ) 



J := 



glog2 



+ 1, (5.4) 



where cq G [1,&o] is any real, and < A < 1, bo > 1 are arbitrary constants. Hence 
J < log and W 1 ^ - Nj x W x / k . We further impose the condition that Hw is 
constant on each subinterval ]Nj,Nj + i], say equal to Hw,j , for j = 0, 1, . . . , J — 1 . By 
simple calculations, 

Nj+i - Nj x 2~ jq W 1/k , (5.5) 



and 

d 2 



dv? 



{-hiyv - « fc ) 1/fc ) = (k - i)hwu k - 2 (w - u k y/ k - 2 x h2^/ k -^ w~ l ' k 



for u G [Nj, Nj + i] . Thus we may apply Lemma 2, part (A), to each of the subintervals 
[Nj, N j+1 ] , with the choice of parameters X = 2~™W 1 / k , Y = h' 1 2-J«( 1 /*-2) w x / k . 
This yields 

E -yh,H Wli E e(-h{W-n k f/ k ) = 

0<h<H W j Nj<n<N j + 1 

W l/(2k) ^ 2^h ( M y/2-i 



— V i h H h y* e (-w 1 ' k (h q + + l ) 

+ 0((\ogW) 2 ) , 



^ ' ()• ;/;< H u 



(5.6) 

since, by the definitions in Lemma 1, jh,H w ^ h 1 . Using the real and imaginary 
parts and taking into account the appropriate portion of (5.2), we readily infer the 
following result. 

Theorem 2. Let W be a large parameter, k > 2 and < A < 1 fixed real numbers, 
Hwj > 1 integers, j = 0, 1, . . . , J — 1 , as explained above, with maxj log Hw,j <^ log W , 
and J given by (5.5) . For a and f3 as defined, and 7 denoting either a or f3 , put 

. , Wl/(2fe) 2 j + 1 h 

Vk-l n ^ . £=20 h 



0<h<H w ^ 

X 



e(w 1/k (h q + £ q ) 1/q - ^) . 
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Then the remainder term Ak(W) defined in (3.3) for the number of lattice points in a 
Lame disc satisfies the estimate 



a»w - 5> (eS;W. w ) 

3=0 

j-i J-i 

< 



E * (Eif WJ ) + E g^±^ + + (log wf) , 

i=o j=o Uw >i 

with an appropriate constant C > . 



Remark. This result may be called a "truncated Hardy's identity for Lame discs". It 
may well be compared with its classic counterpart for the circle, namely 

E r(n)-nX= -X 1 ^ ]T ^ cos(27rV^X-37r/4)+ 0(X l ' 2+£ Y-^ 2 )+0(Y e ) , 

0<n<X 71 l<n<Y U 

where X,Y are any large reals: see I vie [10, f. (1.9)]. Obviously, the formula in 
Theorem 2 is much more complicated. Apart from the technical smoothing factors 
a and (3, this is due to the fact that a Lame curve has points of curvature zero, and 
therefore unavoidable. 



6. Estimation of the multiple exponential sum. Going back to (4.7), it remains 
to bound 

E A fc [{x mk - |m| mfc ) 1/m ) • (6.1) 

|m | <x 

Let U be a further parameter to be chosen later in terms of x. For \ni\ < U, 
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we use (5.1), to get a total error of O ( U (\ogx)~&3w j . The remaining range 

U < n\ < x will be divided into dyadic subintervals in various ways: First, we apply 
Theorem 2, with A = 0.47, to approximate A k {[x mk - |ni| mfe ) 1 / m ). This gives an 
overall error term of 0(x 1Ar ). Next, assuming w.l.o.g. that jj is a power of 2, we 
define U r : =]2 r ~ 1 L r , 2 r U] , r = 1, . . . , R x , where 2 Rl U = f , hence R x < logx. Further, 
V r :=]x(l - 2" r ),x(l - 2" r " 1 )], r = l,...,R 2 , with R 2 such that 2~ R2 ~ 1 x < 10, say. 
Obviously, also R 2 <C logx, and the remaining range gives only a small error of 0(x) . 
After these preparations, and appealing to Prop. 2, it will obviously suffice to bound all 
sums 



EV-(7) 
2^H,W,k,j 



where X is any one of the intervals U r or V r , and < j < J — 1 . Summation over 
r and j will be postponed to the end of the proof. In order to interchange the order 
of summation here, we have to ascertain that H\v,j \ W -( xm k_ nm ky/ m and J do not 
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depend on n\ G X. This is easy for H: Call the corresponding values H(x,j), and 
cf. formulae (6.8) and (6.13) below. It is a bit more delicate for J: Let us observe that 
there exists a constant bo > 1 such that 

max{(x mi - u "»*)i/(mfc) : uel} 
mm{(x mk - u ™fc)i/(™*0 : «gI} < 

for X any one of the intervals U r or V r . Recalling (5.5) and the fact that now 
W = (x mk — n ™ k ') 1 /' m ^ we ma y choose, for every single ri\ G X, 

^mfc j^mfc\l/(mfe) 

CO ~ 



min{(x mA: - W "ifc)i/(mfc) . u e jj • 
Then Co < bo throughout, and J is the same for all ri\ G X. Therefore, 

L>H lXij) ,W,k,j 



VK=(x mfe -n™ fe ) 1 /" 1 



v fe - 1 o<h<H (IiJ - ) £=* h 
(x mk - nT k ) 1/{2mk) e [{x mk - nX lk ) 1/{mk \h q + l q ) 1 ' q } . 



x (6.2) 



x 

ni6X 



Let S(X) denote the complex conjugate of the innermost sum here, and put N = 
(h q + t q ) x l q x £ for short. We transform this sum by Lemma 2, part (B), with 

G(u) = (x mk _ u "»fc)l/(2mfc) j = _jV(x mfc - u "»fc)l/(mfc) _ Apart from the range 

of summation and the weight factors G(n) , this sum is quite similar to that in (5.6), 
with (h, W, k) replaced by (N ', x mk , mk) . Thus Lemma 2, (B), gives 



Y>{Z) = ^UH Y £ qi,2 ~ 1 (N qi + £ qi ) 1/qi ~ s/2 e (-x(N qi +£ qi ) 1/qi + ±) 
Vmk-l ei ^ {x) V J 

+ error terms , 



(6.3) 

with 

mk 

^ 1 mk — 1 

for short. To evaluate the error terms in (6.3), we have to distinguish if X stands for U r 
or V r . Let us call the error terms arising finally in (6.2) S(U r ,j), resp., S(V r ,j). For 
«GW n it follows that G(u) x x 1 ' 2 , and F"(u) x £a: 1 - mfc (2 r t/) mfc - 2 , hence 

G( - U) X £-1/2^/2(^^1-^/2 for ^ £ 

Thus the corresponding error term in (6.3) reads 

O (^- 1 / 2 x mk / 2 {2 r U) 1 - mk/2 ^ +0 (V^logz) . 
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Using this in (6.2) and summing over t and h, we obtain 

S(U r ,j) = O ( x mk ' 2 (2 r U) 1 - mk ' 2 2-M 2 \ogx)+o{Hll A 2 r3) 2-^- 1 ^ 2 x 1 ' 2 logx) , 

(6.4) 

recalling from Lemma 1 that "fh h h is bounded. Similarly, for u G V r , we have 
G(u) x (2-^)V(2mfc) a ,i/2 ; and ^'(w) x£(2- r ) 1/(mfc) " 2 x- 1 , thus 



Therefore, the error terms in (6.3) now read 

O (2-^" 1 / 2 x) +0 ((2" r ) 1/(2mfc) x 1 / 2 logx) . 
Using this in (6.2) and summing over £ and /i, we arrive at 

S(V r ,j) =0 {2~ r 2~^' 2 xlogX^j +Q ^(2- r ) 1 /(2-fc) 2 -j( g -l)/2 (jff - ( ^^ x) l/2 log;z ,j _ 

(6.5) 

It remains to deal with the main term on the right hand side of (6.3), let us call it 
£*(Z). 

Case 1. X = U r . For all v G F'{U r ), we infer that v x i{2 r U) mk ~ x x x ~ mk < £. 
Hence, putting 

E:(X):= £ e (-x(iVHfr) 1/?1 ), 

£i€F'(X), £i<w 

summation by parts readily gives 

£*(W r ) < ^V2-9i/2 x ( £ (2 7 't//x) mfc - 1 )« 1 / 2 - 1 ma_x |£*(W r )| . (6.6) 

We shall estimate E*(W r ) by van der Corput's "fourth derivative test": Writing <f>(w) = 
—x(N qi + w qi ) 1 / qi , straightforward computation gives 

^(w) = - ( qi - 1) x N qi w qi - 4 (N qi + w qi ) 1/qi - 4 x 

x (N 2q i (2 - 91 )(3 - 9i ) + (Mi;) 91 (7 - 4 91 )(1 + qi ) + w 2q ^(l + 9l )(2 + . 

For 1 < 9i < | , which is equivalent to mk > | , the last factor has all nonnegative 
coefficients. For w G F'(U r ) , it thus readily follows that <f>W(w) x i~ 3 x (2 r U/x) 4 ~ 3mk . 
Further, as explained above, length (F'(U r )) <C l{2 r U / x) mk ~ l . In what follows, it will 
be convenient to set 2 r U/x = 2 _s , with 1 < s = Ri — r < R\ — 1. According to [13, 
p. 34, Theorem 2.6], uniformly in v 6W r , 

E*(W r ) < £(2" s ) mfc " 1 (£- 3 x(2- s ) 4 - 3mk ) 1/14 

+ (£(2- s ) mk - 1 ) S/4 (£- 3 x(2- s ) 4 - 3mk )~ 1/14 

^ £11/14 x l/14 |"2-s-j(llmfc-10)/14 _|_ ^27/28 ^-1/14 (2 _s ) ( 2 7mfc-29)/28 
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Using this in (6.6), we conclude that 

£*(W r ) < £ 2/7 x 15/14 (2" s ) 2(mA:+1 )/ 7 + £ 13 / 28 x 13 / 14 (2- s )( 13mfc - 1 )/ 28 
We use this result in (6.2), and sum over I and h, to arrive at 

2^H (UrJ) ,W,k,j 

^ ^2J^(4fe-ll)/(14fc-14) 2-15/14 (2- ss ) 2 ( mk + 1 )/ 7 H 11 / 14 

+ (2 i ) (13fc-27)/(28fc-28) x 13/14 (2 - a) (13mfc-l)/28^/28 + 5( ^ 

where S(U r ,j) has been bounded in (6.4). We balance the first term at the right hand 
side against the error term 0{2~i q 2~ s x 2 H^jj ^) which arises from the last sum in 

Theorem 2. (Recall that U r is of length x 2 r U and ]Nj,N j+1 ] is of length x 2~ jq x: 
cf. (5.5). This gives 



H {u r ,j) 



^ 2 -jyi8fe-ll)/(25fc-25) ^ 2 - s )( 1 0-4mfc)/25 x 13/25j +1 (ggj 



Using this in (6.7) yields 



Tl\ (zUr 



W=(x mk -n™ k ) 1 / r 



;, 7fc + ll 4mfc + 15 w/0 r . ■. 23fc+S4 s 31mfe+35 i^o/mn 

^ (2-J)25( fc -i) (2~ s ) 25 x 37/25 + (2- J )iooc=-D (2~ s ) mo x 143/100 

+ 2-JV2 ( 2 r Uy-mk/2 x mk/2 {Qgx + ^-^^T) (2 r U/x) X 19/25 log X . 

Summing finally over j and r, resp., s, we arrive at**) 
J] A fc ((z mfc - |ni| mfc ) 1/m ) 

U<n x <\x (6.9) 
« x 3 7/25 + jjl-mk/2 x mk/2 bg ^ + t/ -(2m fc -5)/25 x (2m fc +14)/25 bg ^ 

It remains to balance the middle term on the right hand side here against the error term 
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O U x 208 (log x) 8320 j encountered at the beginning of this section. This leads to 



104mfc-131 



U X X 104mfe (log X) 4160 mk 



(*) 

y ' Concerning the last term, we write up this conclusion for 5 — 2mk < . In the contrary case, we 
obtain a term x 19 / 25+£ instead which is negligible. 
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With that choice of U , we finally infer from (6.9) that 
J2 A fc ((z mfc - |ni| mfc ) 1/m ) 

(6.10) 



\m\<2 x 



o 7 / 9 c 339m_fc-262 , 18627mfc - 20614 2238rrafc-655 , _ 

X ' + X 208mfe (logx) 8320mfc + £ 2600mfe T£ 

07/9E 339mfc-262 18627_mfc 

<T ' +j: 208'"'= (lOgx) 8320 



It remains to verify that H(jj r j) , as defined by (6.8), is > 1 throughout. In fact, by 
(5-5), 

H (Urtj) > ( 2 -^)( 18fc - 11 )/( 25fc - 25 )a: 13 / 25 > ^M" - 53 ^" > x o.i384 ^ 
if mk > | . But the case that | < m/c < | brings in only a factor 

» (2- s )Ado-f) »x" 2 / 75 . 

Case 2. X = V r . By a simple calculation, for all v G F'(V r ) , we have u x £2 r / qi . 
Hence, with a look back to (6.3), and S*(V r ) as before, summation by parts yields 

£*(V r ) < 2~ r £- 1 / 2 x max|E;(V r )| . (6.11) 

We estimate again S*(V r ) by the "fourth derivative test". Now 

(4) (w) x f? 1 w" 91 - 3 x x 2" r ( 4 -3/( mfc )) £" 3 x 
for w G V r , hence, again by [13, p. 34, Theorem 6], uniformly in v G V r , 

+ (2^ ^) 3/4 (2-^ 4 - 3 /(^))r 3 x)" 1/14 

^ 2T-(i0mfc-n)/(i4mfc) £ 11 / 14 x 1 / 14 + 2 r ( 29mfc_27 )/( 28mfc ) £ 27 / 28 £ -1 / 14 
Combining this with (6.11), we get 

E* (V r ) « 2~ rJ ^ 11 f/ 7 x 15 / 14 + 2"^^ £ 13 / 28 x 13 / 14 . 
Using this result in (6.2), and summing over t and h, we obtain 

^H (Vrii) ,W,fc,i 

« 2^1^ ^ 4 x i5/i4 ( 6 - 12 ) 
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where S(V r ,j) has been bounded in (6.5). The error term coming from the last sum in 
Theorem 2 now reads 0(2"^2" r x 2 H^ r j) ) , since V r is of length x 2 r x . We balance 
this against the first term on the right hand side of (6.12) - except for the powers of 
2 r , with respect to which we simply compensate the factor 2 r ^smk in the next-to-last 
term of (6.12). This gives 



H, 



(Vr,j) - 



■ 18fc-ll io/nc-1 

2 J25(^tt 2 _r ^^ x 13/25 +1. (6.13) 



With this choice of -ff(v r ,j) > we m fer from (6.12) that 



2^H (VrJ) ,W,k,j 



ni £V r 



+ 0(2-^2-^ 2 ^ 5j) ) 



« ( 2 -i)^y ((2-)^T + (2^) 17/54 ) x 37 / 25 + (2-^)^) * 143 / 100 ((U4) 
+ 2-^/ 2 2- r xlogx + (2- J )^Sr2- r / 54 x 19/25 logx. 
Obviously Hr Vr ^ > 1 throughout: As before after (6.10), 

%,•) » x°- 1384 2^^ » ^.1384-1/27 ^1/10 

Since on the right hand side of (6.14) all exponents of 2 are negative, summation over 
r and j yields 

J2 A fc ((x mfc - Kr fc ) 1/m ) = o(x 37 / 25 ) . 

7jX<|ni |<x 

Together with (6.10) and (4.7), this completes the proof of Theorem 1. 
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